Fracture strength of disordered media: Universality, interactions and tail asymptotics 
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We study the asymptotic properties of fracture strength distributions of disordered elastic media 
by a combination of renormalization group, extreme value theory, and numerical simulation. We 
investigate the validity of the 'weakest-link hypothesis' in the presence of realistic long-ranged inter- 
actions in the random fuse model. Numerical simulations indicate that the fracture strength is well 
described by the Duxbury-Leath-Beale (DLB) distribution which is shown to flow asymptotically 
to the Gumbel distribution. We explore the relation between the extreme value distributions and 
the DLB type asymptotic distributions, and show that the universal extreme value forms may not 
be appropriate to describe the non-universal low-strength tail. 
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It has been known for centuries that larger bodies have 
lower fracture strength. The traditional explanation of 
this size effect is the 'weakest link' hypothesis: the sample 
is envisaged as a set of non-interacting sub- volumes with 
different failure thresholds, and its strength is determined 
by the failure of the weakest region. If the sub-volume 
threshold distribution has a power law tail near zero then 
the strength distribution can be shown to converge to the 
universal Weibull distribution for large sample sizes llj , 
an early application of extreme value theory (EVT) [2j. 

Often failure occurs due to the presence and growth of 
micro-cracks whose long-range interactions call the no- 
tion of independent sub-volumes into question. There 
have been two broad approaches to address such interac- 
tions: fiber bundle models and fracture network models 
Fiber bundles transfer load by various rules as indi- 
vidual fibers fail; in some particular cases exact asymp- 
totic results for the failure distribution have been derived 
[3], and do not explicitly fall into any of the extreme 
value statistics universal forms. Fracture network mod- 
els consider networks of elastic elements with realistic 
long-range interactions and disorder. A particularly sim- 
le approach is based on the random fuse model (RFM) 
[B| , where one approximates continuum elasticity with 
a discretized scalar representation. It has been suggested 
that in the weak disorder limit, fracture would be ruled 
by the longest micro-crack present in the system 
By using critical droplet theory type arguments, one can 
show that an exponential distribution of micro-cracks 
leads to the DLB distribution of failure strengths [7], 
which again does not explicitly have an extreme value 
form. 

These studies raise three important questions. First, 
what is the importance of elastic interactions in deter- 



mining the strength distributions, and does the weak- 
est link hypothesis hold in presence of such interac- 
tions? Second, what is the relation between the DLB 
type asymptotic strength distributions and the universal 
forms predicted by EVT? Third, how should one best 
extrapolate from measured strength distributions to pre- 
dict the probability of rare catastrophic events? We use 
renormalization group (RG) ideas, EVT, and simulations 
of the two dimensional RFM to explore these questions. 
We conclude that (i) the weakest link hypothesis is valid 
for large samples even in the presence of long-ranged elas- 
tic interactions, (ii) the asymptotic forms of the strength 
distribution for these interacting models is compatible 
in disguise with EVT, but of the Gumbel form rather 
than the Weibull form, and, (iii) the use of extreme value 
distributions to estimate the probability of rare events, 
though common in the experimental literature, is not al- 
ways justified theoretically. DLB type asymptotic distri- 
butions (or those derived by Phoenix jij) which depend 
on the details of the material are necessary to safely ex- 
trapolate deep into the tails of the failure distribution. 

The RG and the EVT present two equivalent, yet con- 
trasting, approaches to the study of the universal aspects 
of extreme value distributions in general |lo| . and frac- 
ture strengths in particular. The natural framework to 
investigate the role of interactions and the corrections to 
scaling that emerge as the system size is changed is pro- 
vided by the RG theory. In contrast, the EVT facilitates 
the study of domains of attraction and convergence is- 
sues. The non-universal, yet important, behavior of the 
low reliability tail of the distribution is not described ad- 
equately by either the RG or the EVT. To study such 
non-universal features one needs to develop DLB type 
asymptotic theories. 
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Typically, a RG transformation proceeds in two steps: 
in the first step the system is coarse-grained by eliminat- 
ing short length-scale degrees of freedom, and then the 
resulting system is rescaled. The RG coarse-graining for 
fracture is equivalent to the weakest link hypothesis: a 
system of size L in d = 2 dimensions survives at a stress a 
if its 4 (= 2 d ) sub-systems of size L/2 survive at the same 
stress. This coarse-graining leads to the following recur- 
sion relation for Sl(ct) — the probability that a system 
of size L does not fail under a stress a: 

S L {o)=[S L/2 {a)]\ (1) 

The second step of the RG transformation is to rescale 
the stress suitably and look for a fixed point distribution 
S* that is invariant under RG 

S*(a) =TZ[S*(a)] = [S*(aa + b)} 4 . (2) 

Instead of applying Eq.Q]iteratively like the RG, the EVT 
formulation consider the large length-scale limit directly 

S*(a) = lim [S Lo (Al<J + B^/^ , (3) 

where Lq is a characteristic length-scale. The func- 
tional equations [2] [3] are known to have only three so- 
lutions: the Gumbel, the Weibull, and the Frechet dis- 
tributions. Of these, only the Gumbel (S*(a) — A(er) = 
cxp[— e a ], a £ 3?, a = 1, b = log 4) and the Weibull 
(S*(a) = y a (a) = e-° a , a, a > 0, a = A ( - 1 / a \ b = 0) 
distributions are relevant for fracture. The large length 
norming constants, Al, Bl, satisfy the following asymp- 
totic relations A 2 l/Al — > l/a, \B 2 l — Bl\/Al — > b/a. 

To test the validity of the weakest link hypothesis 
(Eq. [1} in presence of long-range elastic interactions, we 
perform large scale simulations of the RFM 0, [B[ , consid- 
ering a tilted square lattice (diamond lattice) with Lx L 
bonds of unit conductance. Initially we remove a fraction 
1 — p of the fuses at random, where p is varied between 
1 — p = 0.05 and 1 — p = 0.35 (the percolation thresh- 
old for this model is at p = 1/2). Periodic boundary 
conditions are imposed in the horizontal direction and a 
constant voltage difference, V, is applied between the top 
and the bottom of lattice system bus bars. The Kirchhoff 
equations are solved to determine the current distribution 
on the lattice. A fuse breaks irreversibly whenever the 
local current exceed a threshold that we set to one. Each 
time a fuse is broken, we re-calculate the currents in the 
lattice and find the next fuse to break. The process is 
repeated until the system is disconnected. In the present 
simulations, we have considered system sizes from L = 16 
to L — 1024 and various values of p. To explore the 
low strength tail which is beyond the accessible range of 
most experiments, we typically average our results over 
10 5 realizations of the initial disorder. The fuse model is 
equivalent to a scalar elastic problem. Using this equiv- 
alence, the strain is defined as e = V/L and the stress is 
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FIG. 1: Testing the weakest link hypothesis. Comparing the 
survival probability Sl(ct) for a L x L network (solid lines) 
with that predicted by the weakest link hypothesis, Sl/2 ( ct ) 4 , 
(dotted lines) for 1 — p = 0.10. Note the excellent agreement 
even for moderate system sizes. 

given by a = I/L, where / is the current flowing in the 
lattice. The fracture strength is defined as the maximum 
value of <7 during the simulation. 

The RG coarse-graining step (Eq. []} produces a nat- 
ural test for the weakest link hypothesis. In Fig. Q] we 
report the survival probability Sl(&) for different sys- 
tem sizes L, compared with those for systems of size L/2, 
rescaled according to Eq. Q] The agreement between the 
two distributions is almost perfect for L/2 > 32, indi- 
cating that Eq. [TJis satisfied asymptotically. Corrections 
to scaling due to the effect of distant micro-cracks are 
expected to decay as l/£ 2 , as can be shown by a direct 
calculation, but are too small for us to detect in simula- 
tions (Fig. [1]). We also tested wide rectangular systems 
with L x = 2L y , finding larger corrections, scaling roughly 
as 1/L, which are still irrelevant in the large system size 
limit. 

Duxbury et al. related the survival distribution to the 
distribution of micro-crack widths w p}. At the begin- 
ning of the simulation the 'per-site' probability distribu- 
tion of a crack of width w is P(w < w') = 1 — e~ w / w °, 
where wq ~ — 1/ log 2(1— p) |Tl|. Hence, the distribution 
of the longest crack, w m , in a lattice with L 2 sites is given 
by 

P(w m <w')= (l - e -»7V) . (4) 

The stress at the tip of a crack of width w is asymptotic 
to aK^/w, where a is the applied far-field stress, and K 
is a lattice dependent constant. A sample survives until 
the largest crack becomes unstable when its tip stress 
reaches a threshold ath = cK^/w. Therefore, we have 

S L (a)^ (l-e- ( ^ )2 ) L ^D L (a), (5) 
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FIG. 2: Crack width distributions at peak load, 1 — p = 0.10. 
The initial distribution of micro-crack widths (N(w) is the 
number of clusters of width w), is exponential with (dot- 
ted line, slope « — log 2(1 — p)). As the system is loaded, 
a few bonds break before catastrophic failure; these bonds 
usually connect smaller clusters, producing extra cracks at 
large widths. The resulting crack width distribution at the 
peak load exhibits a size-dependent crossover to a different 
exponential slope. Solid lines represent fits to an exponential. 




where cr = (Tth/K y/wo and D L (a) = exp[-L 2 e~( <T °/ <T ) 2 ] 
is the DLB distribution. To apply the above derivation 
to the failure stress, we first check the distribution of 
micro-crack lengths at peak load. As shown in Fig. [2j 
the distribution is exponential, but due to damage accu- 
mulation, the slope of the tail changes with respect to the 
initial distribution. This appears to be due to bridging 
events in which two neighboring cracks join, leading to 
a modification of Eq. [5] as discussed in Ref. [7j. Thus, 
damage accumulation, though very small, is relevant be- 
cause it changes the exponent of the micro-crack distri- 
bution. The exponential form of the crack length distri- 
bution tail, however, suggests that the DLB form should 
still be valid, as demonstrated in Fig. [3] In particular, 
the average failure stress scales as (a) = ctq/ \J\og(L 2 ) 
(Fig. [3k) and the distributions for different L all collapse 
into a straight line when plotted in terms of rescaled co- 
ordinates (Fig. [3Jd). 

Our arguments thus far are seemingly paradoxical. On 
the one hand we have argued on very general grounds 
that the distribution of failure strengths must be cither 
Gumbel or Weibull, while on the other hand we have 
checked that the failure distribution for fuse-networks is 
of the rather different form proposed by Duxbury et al. 
How can this 'paradox' be resolved? It is easy to check 
that the DLB distribution, when rescaled and centered 
properly, yields a Gumbel distribution, i.e., 
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FIG. 3: Testing the DLB distribution of failure stresses, (a) 
The average failure stress as a function of system size L at 
various bond fractions p (symbols) can be fit well by the DLB 
form (solid lines), except close to the percolation threshold 
(1 — p > 0.3). (b) A collapse of the strength distribution for 
different system sizes at 1 — p = 0.1, such that the DLB form 
would collapse onto a straight line. 



as can be demonstrated by a straightforward calculation 
using A L = a /(2(log(L 2 )) 3 / 2 ) and B L = a Q /^i(l?). 
The above result is striking because fracture distribu- 
tions are usually assumed to not be of the Gumbel form, 
since fracture must happen at positive stress, while the 
Gumbel distribution has support for negative arguments 
as well. This is akin to arguing that the normal distribu- 
tion is not valid for test scores since scores must always 
be positive. Nonetheless, it brings us to the issue of con- 
vergence and validity of extreme value distributions as 
opposed to DLB type distributions. 

The extreme value distributions, S*(a) (=A(<r) or 
^(cr)), are a uniform approximation to the true sur- 
vival function, Sl(&), for all a in the limit of large L, 
i.e., 



lim D L (A L a + B L ) = A(<r) 



(6) 



lim 

L— >oo 



sup 



S L (a) - S* 



Br 



A, 



= 0. (7) 
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In contrast, DLB type distributions 12], are based on 



material details, and are asymptotically correct in the 
low reliability tail, i.e., 



lim 



o-^o 1 - Sl(o) 



(8) 



Note that the uniform convergence in Eq. [7] does not 
bound the relative error in the low reliability tail, while 
the asymptotic convergence in Eq. [5] does. 

The above discussion hints at an underlying question: 
How to accurately predict the probability of rare small- 
strength events with limited experimental data? The 
standard practice is to measure the failure distribution of 
construction beams or micro-circuit wires, fit to the uni- 
versal Weibull or Gumbel form, and extrapolate. How- 
ever, as we have argued, this approach can be danger- 
ous. The low reliability tail is non-universal, and must be 
modeled by a theory that, like DLB, accounts for micro- 
scopic details. Such theories, analogous to critical droplet 
theory (low temperatures), instantons (low h), and Lif- 
shitz tails (low disorder, deep in the band gap) are by 
construction accurate in the low reliability tail. It is in- 
teresting to observe that usually the RG and the critical 
droplet theory address continuous and abrupt phase tran- 
sitions, respectively, yet here these two approaches both 
apply to fracture. 

The convergence to extreme value distributions can be 
extremely slow. For the RFM, let z be number of stan- 
dard deviations up to which the Gumbel approximation 
is accurate within a relative error of e. By using the 
Edgeworth type expansions for the extreme value distri- 
butions [131, we find 
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07exp[-^exp[-^exp[...]]], T] < 4e 2 
log?? - 2 log [log 77- 2 log [. . .]], r] > 4e 2 , 



where the ellipsis indicate an infinite recursion, and 
r] = -(4/3)log(l - e)log(L 2 ). For an accuracy of 10% 
at one standard deviation a sample volume of L 2 as 10 18 
is required, while at 2 standard deviations the required 
sample volume is about L 2 s» 10 264 . As a compar- 
ison, for the Gaussian approximation to the mean of 
a sample of M(3> 1) random variables (normalized so 
that E[X] = 0, E[X 2 } = 1, E[X 3 } = 7) we get, 
z ~ A 1 / 3 + A~ 1 / 3 + C'(A- 4 / 3 ), where A = 6e%/M/7, thus 
z w 3 for e = 0.1, M = 3000, 7 = 2, where the value 
7 = 2 corresponds to the standard exponential distribu- 
tion. However, the universal extreme value forms are not 
always dangerous for extrapolation. One can show that 
they are valid asymptotic forms, a la Eq.|H if they satisfy 
the condition of tail equivalence [3, P- 102] [lEj: 



lim 



1-SlQt) 
1-S*(a) 



C, < C < 00. 



(9) 



micro-crack length distribution may be due to the tail 
equivalence of the microscopic and the Weibull distribu- 
tions. 

In conclusion, by using a combination of renormaliza- 
tion group, extreme value theory, and numerical simula- 
tions we have shown that the failure strength of an elastic 
solid with a random distribution of micro-cracks follows 
the DLB distribution which asymptotically falls into the 
Gumbel universality class. The non-universal low reli- 
ability tail of the strength distribution may not be de- 
scribed by the universal extreme value distributions, and 
thus the common practice of fitting experimental data 
to universal forms and extrapolating in the tails is ques- 
tionable. Theories that account for microscopic mecha- 
nisms of failures, the DLB distribution for instance, are 
required for accurate prediction of low strength failures. 
In our study the emergence of a Gumbel distribution of 
fracture strengths is surprising, and brings into question 
the widespread use of the Weibull distribution for fitting 
experimental data. 
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The success of the classical example of a Weibull distri- 
bution of failure strengths emerging from a power-law 



[1 

[2 
[3; 

K 
[5 

[e; 

[7. 
[8 
[9 

[10 
[11 

[12 



W. Weibull, A statistical theory of the strength of mate- 
rials (Stockholm, 1939). 

E. J. Gumbel, Statistics of extremes (Columbia Univer- 
sity Press, New York, 1958). 

M. J. Alava, P. K. V. V. Nukala, and S. Zapperi, Ad- 
vances in Physics 55, 349 (2006). 

S. L. Phoenix and I. J. Beyerlein, Phys. Rev. E 62, 1622 
(2000). 

L. de Arcangelis, S. Redner, and H. J. Herrmann, J. 
Physique Lett. 46, 585 (1985). 

A. M. Freudenthal, in Fracture, edited by H. Liebowitz 
(Academic Press, New York, 1968), p. 591. 

P. M. Duxbury, P. L. Leath, and P. D. Beale, Phys. Rev. 
B 36, 367 (1987). 

P. D. Beale and P. M. Duxbury, Phys. Rev. B 37, 2785 
(1988). 

B. K. Chakrabarti and L. G. Benguigui, Statistical 
physics of fracture and breakdown in disordered systems 
(Oxford University Press, New York, 1997). 

G. Gyorgyi, N. R. Moloney, K. Ozogany, Z. Racz, and 
M. Droz, Phys. Rev. E 81, 041135 (2010). 
The factor 2 in wo ~ — 1/ log 2(1 — p) is due to the fact 
that, on the diamond lattice, there are 2 W crack 'back- 
bones' of width w with a given site as their left end. 
The DLB distribution is asymptotically exact in the limit 



■5 



of 1 — p — > 0. processes (Springer Verlag, New York, 2007). 

[13] L. de Haan and S. I. Resnick, The Annals of Probability [15] C. W. Anderson, Journal of the Royal Statistical Society. 

24, 97 (1996). Series B (Methodological) 40, 197 (1978). 

[14] S. I. Resnick, Extreme values, regular variation and point 



